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Abstract 

Box consistency has been observed to yield exponentially better per- 
formance than chaotic constraint propagation in the interval constraint 
system obtained by decomposing the original expression into primitive 
constraints. The claim was made that the improvement is due to avoiding 
decomposition. In this paper we argue that the improvement is due to 
replacing chaotic iteration by a more structured alternative. 

To this end we distinguish the existing notion of box consistency from 
relational box consistency. We show that from a computational point 
of view it is important to maintain the functional structure in constraint 
systems that are associated with a system of equations. So far, it has only 
been considered computationally important that constraint propagation 
be fair. With the additional structure of functional constraint systems, one 
can define and implement computationally effective, structured, truncated 
constraint propagations. The existing algorithm for box consistency is one 
such. Our results suggest that there are others worth investigating. 

1 Introduction 

Systems of nonlinear equations where the unknowns are reals arise in specialized 
applications such as the study of chemical equilibrium and in robot kinematics. 
A general class of applications of systems of nonlinear equations arises when 
optimizing a function of n variables with multiple local minima. A common 
optimization method is to set the n partial derivatives to zero, and solve the 
resulting set of n equations, which are in general nonlinear. Thus, nonlinear 
equations occur widely in mathematical modeling. 

Until recently, only Newton's method was available for solving such a system. 
This method is good at refining sufficiently good estimates of a solution. Trying 
to use it otherwise is a hit-and-miss affair. The situation was greatly improved 
with the advent of interval arithmetic ||, . 
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A remarkable subsequent development was BNR Prolog Q (later referred 
to as "CLP(Intervals)"), which introduced what came to be called interval con- 
straints. This method can be regarded as an adaptation of the CHIP system Q . 
This constraint processing system associates each variable with a finite set of 
possible values. Instead, in BNR Prolog, the set of possible values is an interval 
of reals. BNR Prolog adopted from CHIP an instance of a constraint propa- 
gation algorithm that turned out to be an instance of Apt's Generic Chaotic 
Iteration algorithm ("GCI" in the sequel) 0]. 

To solve nonlinear systems with CLP (Intervals), one has the advantage of not 
needing derivatives. The equations are decomposed into primitive constraints. 
The resulting constraint system is then subjected to a constraint propagation 
algorithm (the pruning step). The resulting box is split (the branching step), 
whereupon the same is done recursively to the results of the split. By ensuring 
that pruning preserves completeness (that is, does not remove any solutions), 
one ensures that the solving algorithm generates a sequence of boxes that contain 
all solutions. 

Benhamou, McAUester, and Van Hentenryck |^ showed that this simplicity 
comes at a cost: on the Broyden Banded Function, a widely used benchmark, 
exponentially increasing computation time is needed when one only uses con- 
straint propagation for pruning. This is not surprising because of the presence 
of branching. 

By contrast |2j describes Newton, an algorithm that does not require branch- 
ing in this benchmark and only exhibits linearly increasing computation time. 
This remarkable improvement was based on the novel notions of box consistency 
and of the pseudo-zeros used to characterize maximally box-consistent sets. 

Newton avoided branching on this particular benchmark because box-consistency 
achieved stronger pruning than the constraint propagation used in CLP (Intervals). 
However, box consistency was a step back in the sense of only using functional 
interval arithmetic, which contracts only the interval for the function value, 
rather than using the propagation of CLP (Intervals), which has the potential of 
contracting all intervals involved in a relation. Logically, the next improvement 
was to be one that combined the advantages of box consistency with those of 
the relational interval arithmetic used in CLP (Intervals). This step was taken 
by Benhamou, Goualard, Granville, and Puget with their HC4 algorithm 

In this paper, we explore other ideas for improving the use of relational 
interval arithmetic. We describe the use oi probing an interval constraint system 
to improve the bounds obtained by a single constraint propagation. We show 
that this leads to a relational form of box consistency, which is stronger than 
the original notion of box consistency, which we call functional box consistency. 
We show that functional and relational box consistency are but two extremes 
of a spectrum defined by ways of structuring and truncating the iteration in 
constraint propagation. One of these ways can be regarded as the simulation of 
the evaluation of an expression in interval arithmetic. 

Disclaimer Many basic definitions and results need to be covered here. In 
most cases, no attempt at attribution will be made. In the interest of mutual 
compatibility, some definitions are modified. As a result attribution might not 
be welcomed, yet no novelty is involved. Possible novelties are simulation of 
interval arithmetic by constraint propagation, relational box consistency and its 
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computation by probing, and the identification of alternatives to full constraint 
propagation based on structured rather than chaotic iteration. 

2 Constraints and equations 
2.1 Constraint systems 

Definition 1 A constraint system has the following attributes. 

(1) A set {Ti, . . . , r„} of sets called types. 

(2) A set {xi, . . . , Xn\ of variables where Xi is of type Ti for i G {1, . . . , n}. 

(3) A set {Ai, . . . T Am} of constraints where Ai is an atomic formula of first- 
order predicate logic, {xi^ . . . ,Xn} is the set of all variables occurring in {Ai, . . . , Am}. 
For i = 1, . . . , TO, di C {1, . . . , ?i} is such that {xj \ j e di} is the set of variables 
occurring in Ai. 

(4) ^ state, which is Di x • • • x D„, where, for i £ {1, . . . , n}, Di QTi. We say 
that, in this state of the constraint system, Di is the domain of Xi . 

(5) An initial state, which is a state. 

Definition 2 The relation associated with a function or operation f : TZ^ TZ 
is {{xq, ...,Xk) I xo = f{xi, . . .,Xk)}, for k = 0,1,... 

We will say that a function or operation is admissible if the contraction op- 
erator of the associated relation is efficiently computable. This means roughly 
that it can be computed without iteration. Admissible operations include addi- 
tion, subtraction, multiplication, division, maximum, absolute value, power for 
all integer exponents, exp, log, and the trigonometric functions. 

Proposition 1 Let r be the relation associated with an admissible f : TZ'' — > TZ. 
Let the intervals associated with xq, xi, . . . , Xk be [—oo, -l-oo], Ii, . . . , Ik, respec- 
tively. Then the intervals resulting from applying the contraction operator of r 
are, respectively, the intervals f'{Ii, . . . , Ik), Ii, . . . , Ik, where f is the canonical 
set extension of f. 

Definition 3 An equation is E — 0, where E is an expression of type real 
containing only real variables. An equation system consists of a set X of real 
variables and a set of equations containing no variables other than those in X . 

In the conventional way, we consider expressions as trees. The leaf nodes 
are constants or variables; the nonleaf nodes are the operation symbols. 

A distinguishing feature of CLP (Intervals) is that it decomposes equations, 
or other composite expressions, into primitive constraints. These primitive con- 
straints are the relational versions of the building blocks of expressions, which 
are admissible functions. 

It is this decomposition, described in the definition below, that has been 
identified in as the cause for the observed slowness of CLP (Intervals). 

Definition 4 The constraint system C associated with an equation system £ 
depends on a one-one correspondence between the non-leaf nodes of the trees in 
£ and a set of variables that is disjoint from the variables in £ , which is defined 
as follows. 
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Each of the variables in £ also occurs in C, where it is called a "primary 
variable". To each non-leaf node of an expression in £ there corresponds a 
variable in C that does not occur in £ and is called ''auxiliary variable". 

The constraints in C are determined as follows. For every non-leaf node n 
(which is an operation f), with children ni, . . . ,nk, of a tree in £, there is an 
atomic formula in C with variables Xq, xi, . . . , Xk. The predicate in the atomic 
formula denotes the relation associated with f. This formula is a functional 
atom in C. The variables xi, . . . ,Xk are the input variables of the atom; xq is 
its output variable. 

In addition there is in C, for every root of a tree in £ with corresponding 
variable v, an atomic formula v = 0. This is a relational atom in C. 

Definition 5 (Floating-point numbers, intervals) A floating-point number 
is any element of F U {—oo, +00}, where F is a finite set of reals that includes 
0. If X is a finite floating-point number, then x~ (x~^ ) is the greatest (small- 
est) floating-point number smaller (greater) than x. In addition, —00^ — —00, 
—00+ = —M, -\-oo~ — M , and +00+ = +00, where M is the greatest finite 
floating-point number. 

A floating-point interval is a closed connected set of reals, where the bounds, 
in so far as they exist, are floating-point numbers. When we write "interval" 
without qualification, we mean floating-point interval. 

An interval that does not properly contain an interval is called canonical. 

A box is a cartesian product of floating-point intervals. 

Thus canonical intervals are non-empty sets of reals. They may have positive 
width and they may have zero width. Examples are [a~,a], [a, a"*"], and [a, a], 
where a is a finite float-point number. For any real, there is a unique smallest 
canonical floating-point interval containing it. 

In this paper we consider interval constraint systems, which are constraint 
systems where the types are all equal to TZ and where the domains are intervals. 

3 Propagation 

Apt's Generic Chaotic Iteration algorithm (GCI) jl] is of an astonishing sim- 
plicity and wide applicability. The elegance of CLP (Intervals), noted in j2|, is 
due in part to the fact that its constraint propagation is an instance of GCI. 

GCI maintains a pool of operators that still need to be applied. The attrac- 
tion of GCI is that it does not specify any order among these applications. In 
this section we consider orders of application that arc computationally effective 
for constraint systems that are associated with equation systems. We formalize 
application order by means of "traces" , as defined below. 

Definition 6 A trace for a given constraint system C with attributes as in Def- 
inition U has the following components: 

(1) An index sequence t, which is an infinite sequence with elements in 
{1, . . . ,to}. 

(2) A sequence-of-atoms of which the i-th element is the atom At- in C, for 
z-0,1,... 

(3) A sequence-of-constraints of which the i-th element is the relation defined 
by the atom At- in C, for ? = 0, 1, . . . 
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(-i) A sequence-of-contraction operators of which the i-th element is the con- 
traction operator Tt- defined by the atom At- in C, for z = 0, 1, . . . 

(5) A sequence-of-boxes U of which the i-th element is the initial box ofC if 
i = and is Tt._-^{Ui-i) if i>Q. 

As we are primarily interested in the sequence of boxes, we think of the 
sequence of contraction operators as "activations" of the corresponding con- 
straints. We think of the elements of t as "selecting" a constraint to be acti- 
vated. 

The following proposition is based on the fact that the contraction operators 
are monotone nonincreasing and idempotent and that there are finitely many 
domains. 

Proposition 2 See 0, ||/. 

For any interval constraint system with box B as initial state we have: 

(1) The sequence of boxes has a limit for every trace. 

(2) The greatest lower bound of these limits is also a limit of a trace. 

(3) All traces in which the index sequence is fair have the same limit. This limit 
equals the greatest common fixpoint of ti, . . . ,Tm that is less than B. 

(4) This fixpoint is uniquely determined by the constraint system. It is computed 
by a suitable instance of GCI. 

Definition 7 A constraint system is failed ( non- failed ) if its fixpoint is empty 
(non-empty). 

A failed constraint system has no solutions. A non- failed constraint system 
may, but need not, have solutions. 

Proposition 3 The fixpoint of a constraint system contains all its solutions. 

Definition 8 A segment of a trace is functional if 

(1) The sequence- of- atoms only contains functional atoms. 

(2) For every atom, any input variable that is an auxiliary variable has 
occurred as output variable earlier in the segment. 

(3) No atom occurs more than once. 

A segment of a trace is inverse functional if every occurrence of a variable as 
an output variable has been preceded by an occurrence of that variable as input 
variable or as variable in an activation of an equality constraint. 

The following proposition shows that a trace of a constraint system can 
simulate the evaluation of an expression. 

Proposition 4 Let C be an interval constraint system with attributes as in Def- 
inition U that is associated with an equation system £ containing an expression 
E. Let box B be the initial state ofC such that all its projections corresponding to 
auxiliary variables are [~oo, +oo\. Let d C {1, . . . , n} be such that {xj | j G d} 
is the set of variables in E . Let Xi be the variable in C that corresponds to the 
root of E. Let T be a functional initial segment of a trace. 

The i-th projection of the last of the sequence-of-boxes of T is the same 
interval as the one obtained when E is evaluated in interval arithmetic with 
TTj{B) as the interval substituted for xj in E, for all j G d. 



5 



Consequences of Proposition |j GCI leaves open the possibility of activat- 
ing an operator that has no effect; that is, when this activation does not result in 
any domain reduction. In any trace containing an operator activation without 
effect, the number of steps to convergence can be decreased by removing it. 

One possible heuristic to avoid vacuous activation is to require a trace to be a 
two-phase iteration: to consist of a repetition of cycles consisting of the following 
two items: (1) a functional segment and (2) an inverse-functional segment. 

It may be conjectured that such a trace is optimal in the sense of having a 
shortest pre-convergence initial segment. It is to be expected that such a trace 
is subject to severely diminishing returns in the sense of activations resulting in 
domain reduction. 

The first functional segment likely has the greatest effect. It corresponds to 
the interval-arithmetical special case of constraint propagation. While one can 
probably construct examples to the contrary, most of the remaining reduction 
is typically effected in the first inverse- functional segment. For some purposes, 
to be discussed below, it is advisable to truncate constraint propagation before 
convergence has occurred. Promising truncations are: after the first functional 
segment or after the first cycle. 

4 Box consistency 

We first describe, adapted to the current setting, the box consistency notion of 
. To distinguish it from the version to be described next, we call it "functional 
box consistency" . 

4.1 Functional box consistency 

Definition 9 (Coordinate-wise functional box consistency operator) 

Given an equation E — with variables Xi, . . . , Xn, and i G {1, . . . , n}. The i-th 
coordinate-wise functional box consistency operator replaces m /i x • ■ • x /„ the 
interval li by its intersection with the smallest interval containing {y £ 7?. | G 
-E'}, where E' is the result of evaluating E in interval arithmetic with /i, . . . , /„ 
substituted for the variables xi,. . . ,Xn, except that the smallest floating-point 
interval containing y is substituted for Xi. 

The i-ih coordinate- wise functional box consistency operator is completeness- 
preserving in the sense that it removes no solution to E = Q. 

Definition 10 (Functional box consistency) A box is functionally box con- 
sistent with respect to an equation if it is a common fixpoint of the n coordinate- 
wise functional box consistency operators associated with that equation. 

This is equivalent to the notion of box consistency introduced in Q . By itself, 
box consistency is not interesting: for example, the empty box is box-consistent 
with respect to {x = 0}, which has plenty of zeros. What is lacking so far is a 
suitable notion of maximality. This was done in |^ by relating box consistency 
with leftmost and rightmost pseudo-zeros. It can also be done by means of 
fixpoint theory. Because each coordinate-wise functional box consistency oper- 
ator is a monotonic one defined on the partially ordered set of boxes with set 
containment as partial order, it has a greatest fixpoint. By applying GCI one 
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obtains the greatest fixpoint common to all n coordinate-wise functional box 
consistency operators associated with a set of equations. 

Definition 11 (Functional box consistency operator) Given an equation 
E = with variables xi, . . . ,Xn- The functional box consistency operator is the 
mapping from a box B to the greatest fixpoint contained in B that is common 
to all n coordinate-wise functional box consistency operators associated with the 
equation. 

Definition 12 (Functional pseudo-zero) Given an expression E with vari- 
ables xi, . . . ,Xn. An i-th functional pseudo-zero of E with respect to the inter- 
vals {Ij I j e {1, . . . , n} \ {i}} is a canonical interval I such that is contained 
in the interval resulting from evaluating E in interval arithmetic with I substi- 
tuted for Xi and Ij substituted for xj, for all j € {1, . . . ,n}\ {i}}. 
An i-th functional pseudo-zero of an equation system with respect to the in- 
tervals {Ij I j e {1, . . . ,n} \ {i}} is a canonical interval that is a functional 
pseudo-zero of every one of its equations. 

Proposition 5 Given an expression E with variables xi, . . . ,Xn. If, for any 
i G {1, . . . ,n}, no i-th functional pseudo-zero exists with respect to {Ij \ j E 
{1, . . . , ^'^en E = has no solution with Xj £ Ij, for allj e {!,..., n}\ 

Proposition 6 The result of applying the functional box consistency operator 
of an equation system with at least one solution to a box B results in a box that 
has as i-th projection the least interval containing the i-th leftmost and rightmost 
functional pseudo-zeros with respect to tti (B), . . . , iTniB), for alii e {1, . . . ,n}. 

4.2 Probing 

Definition 13 Let C be a constraint system. Probing C with a constraint A 

means determining whether the constraint system C , which is the result of 
adding A to C, is failed or non-failed. 

Note that C does not change as the result of probing. 

Proposition 7 (Monotonicity of probing) // a non-failed constraint sys- 
tem C yields failure as a result of probing with x <ui, then it also yields failure 
as a result of probing with x < U2, where ui and U2 are floating-point numbers 
such that U2 < ui . 

A similar fact holds for probing with x > h ot x > I2, with I2 > h. Also non- 
failure can be inferred from non-failure by means of suitably selected probes. 

When probing a constraint system C containing a variable x with x < a 
yields failure, then it has been proved that no solution has an a;-component 
that is less than or equal to a. It does not follow that probing C with x > a 
yields non-failure. It is quite common for C to be non-failed, yet not to have 
any solutions. This is sometimes discovered by probing C with both x < a and 
X > a and finding failure in both cases. While this property of probing is often 
effective, such an x and such an a cannot always be found. 
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As seen above, probing has a logic of its own. Here is an other example. 
It may be that a non-failed constraint system C containing variable x yields 
non- failure on probing with x > a and that it yields non-failure on probing with 
X < b, where a < b, and yet yields failure on probing with x > a A x < b. 
Conversely, non-failure on probing with x > a A x < b implies non-failure on 
probing with a; > a; it also implies non-failure on probing with x < b. 

Probing suggests a relational version of the functional pseudo-zero. As prob- 
ing applies to any constraint system, not just to those derived from equations, 
we call them "pseudo-solutions" . 



4.3 Relational box consistency 

Probing can be used to compute approximations to relational box consistency, 
a criterion similar to functional box consistency, but one that produces closer 
approximations to the set of solutions. Another way in which relational box 
consistency is interesting is that it applies to all interval constraint systems, not 
only to those that are derived from equations. 



We proceed in analogy with Section 4.1. In analogy with the pseudo-zeros 
of we have the following definition. 

Definition 14 ((Canonical) pseudo-solution) Let C be a non-failed inter- 
val constraint system with initial box /i x • • • x /„ and let i S {1, . . . ,n}. Let 
C be an interval constraint system with the same attributes as C, except that 
interval Li is changed to an interval y. If C is non-failed, then y is an i-th 
pseudo-solution of C with respect to {Ij \ j € {1, . . . ,n} \ {i}}. If y is a canon- 
ical interval, then it is an i-th canonical pseudo-solution of C with respect to 
{/, |je{l,...,n}\{*}}. 

Functional pseudo-zeros have been defined, following [|[ ^ as being canonical 
intervals. Here it is useful that pseudo-solutions are not necessarily canonical. 

Definition 15 (Coordinate-wise relational box consistency operator) 

The i-th coordinate- wise relational box consistency operator of an interval con- 
straint system C maps the initial box /i x • ■ • x /„ to one where the i-th pro- 
jection (i G {!,..., n}j has been replaced by the least floating-point interval 
containing the union of the i-th canonical pseudo-solutions of C with respect to 
{J, |je {!,..., n}\W}. 

We introduced Definition |l^ because it figures so prominently in the litera- 
ture, following Q. One could define a relational analog, but it is more useful to 
skip to the analogs of Definition |ri| and Definition |l^. 

Definition 16 (Relational box consistency operator) Given an interval con- 
straint system C. The relational box consistency operator of C maps its initial 
box /i X • ■ • X J„ to the greatest flxpoint contained in it that is common to all n 
coordinate-wise relational box consistency operators associated with C. 



Definition 17 (Leftmost (Rightmost) pseudo-solution) Let C be an in- 
terval constraint system with variables xi, . . . ,Xn and let i G {1, . . . , n}. An i-th 
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pseudo- solution [a,b] of C is an i-th leftmost pseudo-solution of C if C has no 
solution (^1, . . . , with in {x d TZ \ x < y for all y 6 [a, 6]j£]. 

Pseudo-solutions become less interesting the wider they are. For example, 
[—00, -l-oo] is a pseudo-solution, a leftmost, and a rightmost pseudo-solution for 
any interval constraint system. On the other hand, consider an i-th canonical 
pseudo-solution [a, b] such that no z-th canonical pseudo-solution is to the left 
of it. As every real that is the i-th component of a solution is contained in 
an i-th canonical pseudo-solution, [a, b] is also an i-th leftmost pseudo-solution. 
Hence the following proposition. 

Proposition 8 Suppose the relational box consistency operator maps box B 
to B' . Then B' is equal to the smallest box containing the i-th leftmost and 
rightmost canonical pseudo-solutions, for i G {1, . . . , n}. 

The reason for having these different characterizations for the same object 
is that the one in terms of pseudo-zeros is convenient for a relational box- 
consistency algorithm; see section ^. The one in terms of fixpoints is convenient 
for comparing relational box consistency with functional box consistency. 

Proposition 9 Let E be an expression with n variables and let C be the interval 
constraint system associated with {E — 0}. The result of mapping a box B with 
the relational box consistency operator of C is contained in the result of mapping 
B with the functional box consistency operator of {E — 0}. 

Proof. Both mapping results are characterized as common greatest fix- 
points. Hence it is sufficient to show that, with respect to {Ij \ j £ {1, . . . , n} \ 
{«}}, the result of the i-th relational box consistency operator is contained in 
that of the i-th functional box consistency operator. 

This can be ascertained by considering a canonical interval [a, 6] for which 
^ fi{[a,b]). Here fi is a function from intervals to intervals such that fi{I) is 
the result of evaluating E in interval arithmetic with Ij substituted for Xj for 
j G {1, . . . , n} except that I is substituted for Xi. 

Such an interval [a, b] lies outside the result of applying the i-th functional 
box consistency operator. One fair trace of C begins with a functional segment 
followed by activating the equality constraint. It yields the empty interval at 
that point, which ensures that the unique limit of any fair trace is empty. Hence 
[a, b] also lies outside the result of the i-th relational box consistency operator. 

5 Computing Box Consistency 

As only maximal box-consistent boxes are of interest, we compute leftmost and 
rightmost pseudo-zeros. 

5.1 Computing the leftmost functional pseudo-zero 

Given an equation E with variables .ti, . . . , Xn- The purpose of function ZEROl 
described in this section is to compute the i-th functional leftmost pseudo-zero 
of E. It is an adaptation of function LeftNarrow in Q. 

^ It is simpler to say: "with < a", but this does not yield a useful definition because 
is a real and a is a floating-point number. As a result, a can be —00. 
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interval function ZERO 1 (a, 5) { 
if (0^/([a,6])) return 0; 
if {[a,b] is canonical) return [a,b]; 
//[a,b] not canonical, so has midpoint 
m := midpoint of [a, 6]; 
/' := ZEROl(a, m); 
if (/' is empty) return ZEROl(m, b); 
return /'; 



Figure 1: A definition of a function to compute a functional leftmost pseudo- 
zero. 

interval function ZER02(a,6) { 

if ([a, 6] is canonical) return [a, 6]; 

m :~ midpoint of [a, 6]; 

probe C with Xi < m; 

if (result of probing is failure) { 

probe C with m < Xi A Xi < b; 

if (result of probing is failure) return 0; 

return ZER02(m, 6); 

} 

/' ;~ ZER02(a, m); 

if (/' is empty) return zero2(to, 6); 

return /'; 



Figure 2: A definition of a function to compute the i-th canonical leftmost 
pseudo-solution. 



Suppose function ZEROI is first called with arguments ag and bo such that the 
j-th leftmost functional pseudo-zero, if it exists, is contained in [ag, bo]. It returns 
an interval [a', 5']. If this interval is empty, the algorithm has detected that 
no solution occurs in [ao,6o]. Otherwise, [a', 6'] is the i-th functional leftmost 
pseudo-zero. 

See Figure ^.l| for a definition of zeroI. Here /([a, 5]) denotes the result of 
evaluating E with Bj substituted for variable Xj for all j G {1, . . . , 71} \ {i} and 
[a, 6] substituted for Xi. 



5.2 Computing canonical leftmost pseudo-solutions 

The function zero2 is defi ned recursively to compute the relational leftmost 
pseudo-zero; see Figure 5.2 . 

It assumes and maintains the invariant that when zero2 is called with a 
and b as arguments, the interval [a,b\ is an i-th leftmost pseudo-solution. The 
result is when it has been proved that no solution exists in [a,b]. If the result 
is not 0, then it is [00,^*0] such that [00,^*0] ^ [oj and oq < &0j and [ao,6o] is 
the i-th canonical leftmost pseudo-solution. 
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Definition 18 (Functionally truncated probing) Let C be the constraint 
system associated with an equation system. Functionally truncated probing of C 
with a constraint A refers to the result (failure or nonfailure ) of an initial func- 
tional segment of a trace of the constraint system C that results from adding 
A to C. The result is nonfailure if is contained in all intervals of variables 
occurring in the equality constraints in C ; failure otherwise. 

If in ZERO 2 one would replace probing with functionally truncated prob- 
ing, then an algorithm would result that is for practical purposes equivalent 
to ZEROl. That is, functional box consistency can be computed with the same 
efficiency by constraint propagation on a system of atomic constraints, provided 
that propagation is not chaotic, but suitably structured and truncated. 

6 Conclusions 

Following lH, we have treated systems of equations. Functional box consis- 
tency is easy to generalize to systems containing both equalities and inequalities. 
Relational box consistency seems more general because it applies to all interval 
constraint systems, not just to those that are derived from systems of equalities 
and inequalities. 

In 1^ Newton was compared with CLP (Intervals) on the Broyden Banded 
function. CLP (Intervals) was observed to required time exponential in the num- 
ber of variables, whereas Newton required linear time. This is indeed to be ex- 
pected: CLP (Intervals) used for pruning a single application of GCI. Because of 
the weakness of such pruning, the search tree reaches a significant depth. The 
size of such a tree is exponential in the number of dimensions. 

In Q, the observation was made that Newton requires no branching on this 
example. Hence no exponential behaviour is to be expected. 

It is now time to look beyond this particular example to those where even 
with pruning as powerful as in Newton, substantial branching is necessary. 
Newton showed that more effort spent in pruning is rewarded by a reduction in 
branching in such a way that the total computation time is much reduced. Thus 
there is a trade-off between time spent on pruning and time spent on branching: 
at some point, additional effort spent on pruning must stop being productive. 
This might suggest replacing functional by relational box consistency. 

There is a better method. Note that to compute both functional and re- 
lational box consistency, one iterates all the way down to canonical intervals, 
the narrowest that the floating-point hardware allows. This is done to make 
the interval for one variable as narrow as possible. Yet to compute box consis- 
tency one has to do this for all variables repeatedly until no further narrowing 
is possible for any variable. In the beginning, most the intervals for most of 
the variables are still wide. While this is the case, it seems wasteful to iter- 
ate in functions zeroI or zero2 all the way down to canonical intervals: the 
convergence criterion should be adapted to the width of the other intervals. 

In addition to this improvement, which applies both to ZEROl and to zero2, 
there is an improvement that applies to the latter alone. As we described probing 
here, propagation is completed to convergence. As shown in Proposition ^ 
truncating the trace in propagation to the initial functional segment, causes 
ZER02 to compute functional box consistency. By truncating the trace less 
drastically, say, till after the first cycle of a two-phase iteration, one obtains 
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better chance at getting failure in probing, yet avoids the neghgible reductions 
associated with the last phases of propagation to convergence. It seems worth 
investigating how many phases are optimal in this respect. 

7 Related work 

The routine absolve in BNR Prolog uses probing to find narrower intervals 
than a single application of GCI can give. The mechanism was discovered inde- 
pendently by Chen and van Emden , who called it "hypernarrowing" . It uses 
bisection to determine the greatest m such that probing with x < m gives fail- 
ure. In Q "hypernarrowing" was used in optimization. A dramatic decrease in 
the number of function evaluations resulted. missed the connection between 
"hypernarrowing" and box consistency. 

Benhamou et al. 1^ noted the ineffectiveness of the CLP (Intervals) solve 
routine for nonlinear equations. In response they introduced box consistency 
and used it to achieve dramatic improvements over the CLP (Intervals) solve. 
In the version they introduced (here called functional box consistency), they 
implicitly discard constraint propagation, and only use interval arithmetic. 

The HC4 algorithm of || is a propagation algorithm where instead of indi- 
vidual contractions one applies an algorithm called HCRrevise, which similar to 
a two-phase iteration truncated after the first cycle. 
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